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Constraints on Black Hole Spin in a Sample of Broad Iron Line AGN 

Laura W. Brenneman 1 , Christopher S. Reynolds 2 

ABSTRACT 

We present a uniform X-ray spectral analysis of nine type-1 active galactic nuclei 
(AGN) that have been previously found to harbor relativistically broadened iron emis- 
sion lines. We show that the need for relativistic effects in the spectrum is robust even 
when one includes continuum “reflection” from the accretion disk. We then proceed to 
model these relativistic effects in order to constrain the spin of the supermassive black 
holes in these AGN. Our principal assumption, supported by recent simulations of 
geometrically-thin accretion disks, is that no iron line emission (or any associated X- 
ray reflection features) can originate from the disk within the innermost stable circular 
orbit. Under this assumption, which tends to lead to constraints in the form of lower 
limits on the spin parameter, we obtain non-trivial spin constraints on five AGN. The 
spin parameters of these sources range from moderate (a ~ 0.6) to high (a > 0.96). 
Our results allow, for the first time, an observational constraint on the spin distribu- 
tion function of local supermassive black holes. Parameterizing this as a power-law 
in dimensionless spin parameter (/(a) °c |a|^), we present the probability distribution 
for C, implied by our results. Our results suggest 90% and 95% confidence limits of 
C, > —0.09 and £ > —0.3 respectively. 

Subject headings: accretion, accretion disks - black holes - galaxies:nuclei - X- 
rays: spectra 


1. Introduction 

Black holes (BHs) are among the simplest objects in nature, able to be defined completely 
by mass and spin alone. But whereas the mass of a BH is relatively easy to calculate, provided 
material orbiting the BH can be observed, spin has proven much more challenging to constrain. 
The advent of high-resolution X-ray spectroscopy and the promising development of gravitational 
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wave detectors are now making this field of research accessible. This paper represents the first- 
ever survey to quantify the angular momenta of BHs in a sample of AGN using relativistically 
broadened iron lines as spin diagnostics. We anticipate that this work will serve as a stepping stone 
for guiding future spin surveys across a wide range of environments, including AGN of varying 
host galaxy morphology, and supermassive BHs vs. stellar-mass BHs, to name two interesting 
population comparisons. 

A General Relativistic effect, BH spin manifests itself physically by setting the radius of the 
effective inner edge of the accretion disk: as the magnitude of the (prograde) spin increases, frame- 
dragging provides matter in the disk with additional angular momentum with which to resist the 
inward pull of gravity, enabling this material to remain in orbit nearer to the event horizon than it 
would in such close proximity to a non-spinning BH. The inner edge of the disk can therefore be 
used to measure BH spin (Reynolds & Fabian 2008). 

The position of the inner edge of the accretion disk can be constrained through detailed mod- 
eling of the X-ray spectrum of the disk. Two popular modeling schemes include fitting the thermal 
continuum of the disk (Narayan et al. 2007) and fitting the shape of the fluorescent Fe-K line emit- 
ted by the disk in response to irradiation by the X-ray power-law (Brenneman & Reynolds 2006). 
The former method relies on a priori knowledge of the BH mass and distance as well as the disk 
inclination angle, however, and is also only applicable to smaller Galactic BH systems (GBHs) in 
thermally dominant states as opposed to supermassive BHs found in the centers of active galac- 
tic nuclei (AGN). By constrast, the shape of the iron line profile depends upon the gravitational 
redshift and fluid velocity field, both of which depend only upon M/R. The iron line profile is 
therefore mass-independent, making it a useful technique for any BH system (Reynolds & Nowak 
2003), but especially for AGN in which the NH mass is often poorly determined. 

The shape of the iron line profile is determined by several physical parameters: the BH angu- 
lar momentum, the radial extent of the disk, the emissivity profile of the disk, the inclination angle 
of the disk to the observer’s line of sight and the flux in the line itself as compared to the contin- 
uum flux. A proper model for the iron line must also take into account the effects of the Doppler 
shift. Special Relativity (e.g., beaming) and General Relativity (e.g., gravitational redshift). Un- 
fortunately, the two widely available models within the xspec spectral analysis software for these 
types of fluorescent disk lines emitted near a BH do not take all of these physical effects and pa- 
rameters into account. BH spin, in particular, is not a free parameter in either the diskline model 
for non-spinning Schwarzschild BHs (Fabian 1989) or the laor model for maximally-rotating Kerr 
BHs (Laor 1991). Both of these codes have their dimensionless BH spins hardwired at a = 0 and 
a = 0.998, respectively, where a = cJ/GM 2 . If one wishes to constrain BH spin, therefore, it is 
necessary to create a new, fully-relativistic model code which enables the BH spin to be fit as a 
free parameter. 
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We have developed such a model as both a line code (kerrdisk) and a convolution kernel 
(kerrconv) that can be applied to a full spectrum of fluorescent emission lines from the disk. 
These models were introduced and detailed in Brenneman & Reynolds (2006), hereafter known as 
BR06. In that paper we also presented a series of model fits to the 330 ks XMM-Newton/EPIC-pn 
spectrum of the canonical broad iron line Sy-1 AGN MCG-6-30-15, noting that the best fit was 
achieved with a kerrconv model acting on an ionized disk reflection spectrum, where the inner 
edge of the disk was r m i n < 1.62r g . These data have allowed us to place formal constraints on the 
BH spin in MCG-6-30-15 of a > 0.987 (90% confidence level). Because this fit was performed 
under the assumption of a strict truncation of the iron line region at the innermost stable circular 
orbit (ISCO) in the disk, this constraint is somewhat weakened ( a > 0.94) if the iron line emission 
region does in fact bleed inside of the ISCO (Reynolds & Fabian 2008). Given that MGC-6-30-15 
harbors the broadest iron line in an AGN observed to date, this model fit reinforces our expectation 
that this AGN contains a BH of near-maximal spin. Many other, similar Sy-1 galaxies are also 
observed to have broad iron lines, however (Guainazzi et al. 2006; Nandra et al. 2007; Miller 
2007), which begs two important questions: (1) can XMM-Newton observations of broad iron lines 
be used to robustly constrain BH spin, and (2) how is spin distributed between and among different 
types of BH systems? 

To begin addressing these questions, we have collected a sample of eight other Sy-1 AGN 
from the XMM-Newton archive that have previously been observed to have broad iron lines. Us- 
ing the same technique we used to model the spectrum of MCG-6-30-15 by first addressing the 
continuum and absorption and then considering the broad line region, we assess each AGN inde- 
pendently, with the goal of obtaining statistically robust BH spin constraints whenever possible. 

The question of what BH spin distribution we should expect to see among Sy-1 AGN is not 
a simple one. Unlike GBHs, whose spins are thought to be chiefly natal due to the small masses 
and/or short lifetimes of their stellar companions, AGN are old and massive, and are thought to 
have grown from a combination of mergers and accretion over their long histories. By compar- 
ing observations of BHs in local galaxies and AGN to a model of what the BH mass function 
would look like for a given local population whose BHs were grown solely by accretion, Mar- 
coni et al. (2004) find that mergers are not significant in the evolution of local supermassive BHs. 
The authors suggest that mergers might play a more prominent role for SMBHs beyond z = 3. 
Determining the spin of such sources would therefore provide us with a new parameter to use in 
constraining the role and efficiency of accretion in evolving local SMBHs, and would also aid in 
assessing the relative contributions of accretion vs. mergers in higher redshift sources. Similarly, 
Volonteri et al. (2005) have used numerical simulations to predict how massive BHs should form 
and evolve. They find that the coalescence of comparable-mass BHs increases their BH spins, 
while the capture of smaller BH companions in randomly oriented orbits reduces the spin of the 
primary BH. Given the distribution of massive BH binary mass ratios in hierarchical models, binary 
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mergers alone do not lead to a systematic spin-up or spin-down of SMBHs with time. By contrast, 
because of the alignment of a SMBH with the angular momentum of the outer accretion disk, gas 
accretion tends to increase BH spin over time. Overall in the Volonteri et al. models, accretion 
dominates over mergers and efficiently increases BH spin. The spin distribution is heavily skewed 
toward rapidly rotating holes, is already in place at early epochs, and does not change much below 
z = 5. If accretion occurs via a geometrically thin disk, about 70% of all massive BHs are close 
to maximally rotating. Even in the conservative case in which accretion is via a geometrically 
thick disk, about 80% of all massive BHs have a > 0.8 (Volonteri et al. 2005). King, Pringle & 
Hofmann (2008), however, use simulations to contend that black hole growth in AGN proceeds by 
accretion from random directions such that the average BH spin ranges from a = 0.1— 0.3, with 
mergers making little difference. These authors argue that more massive BHs would tend to have 
smaller spins, and that the “rare” examples of massive BHs with large spins would likely result 
from prograde mergers of SMBHs with similar masses and would therefore be most often found 
in giant elliptical host galaxies (King et al. 2008). Given these disparate theoretical predictions 
and the increasing computing power of complex numerical simulations, the time is ripe for culling 
observational evidence to investigate whether the data support such hypotheses. 

In order to begin probing the question of BH spin distribution across the AGN population we 
clearly need to expand our sample. In this paper we restrict our attention to several other nearby 
Sy-1 AGN in an effort to reduce the number of environmental variables that may be involved in 
spin determination (e.g., no GBHs, though we do intend to expand our source list further in the 
future to accommodate such objects and compare the results to those of our AGN sample). Because 
the kerrdisk and kerrconv models rely on the presence of broad iron lines in the X-ray spectrum 
in order to constrain BH spin, we have also restricted our source list to include only those Sy-1 
AGN that have robustly been observed to harbor broad iron lines based on previous studies with 
XMM-Newton and/or Chandra (Nandra et al. 2007; Miller 2007). 

Section 2 discusses the observations of the sources in our sample and the method used to 
reduce the EPIC-pn data and extract spectra. A summary of our modeling technique is given in 
§3. Results and tables comparing the relevant best-fit parameters for each source can be found in 
§4. Implications of these results, conclusions we have drawn and a prognosis for future work in 
this field are discussed in §5. The Appendix contains the detailed spectral fitting results for each 
source. 


2. Observations and Data Reduction 

Our source selection was based on recent work by Nandra et al. (2007) and the review by 
Miller (2007), both of whom consider the robustness of broad iron line detections across different 
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source populations. The Miller review, in particular, collects all the recent X-ray data published on 
Seyfert sources and categorizes these objects by the strength and robustness of their broad iron line 
features. We have used these two studies as a guide for selecting the candidates for our sample. 

Our data have been collected from the XMM-Newton Science Archive, using only publicly 
available observations. XMM-Newton data were preferentially used as opposed to Chandra data 
because of the superior collecting area of XMM-Newton in the 2 — lOkeV band where the broad 
Fe-Ka feature is the most prominent. Though Suzaku is producing very interesting results, there 
are not enough publicly available data collected by that instrument to undertake such a study at 
this time. 

We focus our study exclusively on the EPIC-pn data to avoid cross-calibration issues between 
the pn, MOS and RGS instruments. Though the RGS spectra in particular can provide detailed 
insight into the 0.3 — 2.0keV spectrum and help illuminate the nature of the soft absorption and 
excess seen in so many Sy-1 AGN, the cross-calibration issues unfortunately render it nearly im- 
possible to correlate those results with the continuum and iron line spectrum observed by the 
EPIC-pn instrument. For further discussion of this issue, see Brenneman et al. (2007) and refer- 
ences therein. Because of the higher effective area of the pn, these data are more useful in our 
work. Observation data for each source is shown in Table 1. 

For each source in question we performed reprocessing (when necessary) and data reduction 
with the SAS version 7.0.0 software, including the latest CCF calibration files. We have used 
the SAS epatplot task to compute the fraction of single, double, triple and quadruple events 
as a function of energy and compared these fractions to their nominal values as measured from 
weak source observations. For sources that are affected by pile up, these fractions deviate from 
the nominal values due to the higher probability of wrong pattern classification. No significant 
deviation from the nominal single and double distributions was found for any of our sources, 
indicating that our EPIC-pn observations are not affected by pile up. Source and background 
spectra and light curves were extracted using the xmms elect task from the SAS GUI. Response 
matrices and ancillary response files were generated using the rmf gen and arf gen tasks. These 
were then grouped with the spectral files using the FTOOLS package grppha with a minimum of 25 
counts/bin in order to use x 2 a s a meaningful goodness-of-fit statistic. Spectral analysis was carried 
out as before for MCG-6-30-15 (BR06) with XSPEC version 11.3.2 and our local kerrdisk and 
kerrconv models installed. 1 


'These models are available for download and use in XSPEC from 
http://www.astro.umd.edu/chris/publications/papers/kerrdisk/kerrdiskjnodel.html. 
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3. Spectral Analysis Methodology 

In the following Section we outline the method used for the data reduction and spectral anal- 
ysis of the AGN in our sample. We then compare the results from all the sources in our sample 
and discuss their physical implications for BH spin distribution and evolution among Sy-1 AGN. 
For spectral modeling and results from the individual sources in our sample, see the Appendix. 


3.1. Continuum Fitting 

We have endeavored to make the spectral analysis method as uniform as possible for all the 
AGN examined. There are, however, certain intrinsic physical differences between the sources that 
must be taken into account in order to properly model the continuum and isolate the iron line(s) 
in each case. Though all of our sources are Sy-1 AGN (or otherwise observed in states where 
they manifest broad disk lines), each system is unique in its physical properties. Some have more 
photons than others to work with, either because of higher flux, longer observation time or both. 
Some exhibit evidence for complex, multi-zone warm absorbers (WAs) intrinsic to the central 
engine, whereas some show only cold absorption from neutral hydrogen along our line of sight. A 
soft excess is seen in some sources but not others. And of course, the strength and breadth of the 
Fe-Ka line varies from source to source as well, though all have been chosen because they have 
previously shown evidence for broad iron lines. 

In order to examine the iron line(s) in detail the continuum must first be accurately modeled. 
After excluding the energy ranges relevant to the iron line (~ 4.0 — 8.0 keV) and the mirror edges 
(~ 1.5 — 2.5 keV), we examine the rest of the 0.6 — 10.0 keV spectrum, fitting it with a power- 
law continuum typical of an AGN which is then modified by cold photoabsorption from neutral 
hydrogen. We set the minimum value of jVh equal to the Galactic column density along our line 
of sight to the source 2 and allow it to vary as necessary to accommodate the cold absorption in 
the system. Some sources require a second, unabsorbed power-law component as well to properly 
model the basic continuum shape. We disregard energies below ~ 0.6 keV due to calibration 
uncertainties in this range for the EPIC-pn instrument. 

If significant residual features remain after fitting this absorbed power-law which indicate 
the presence of a soft excess and/or warm absorption, we include these components one by one 
as long as they make a significant difference in the global goodness-of-fit according to the F-test 
(Ax 2 < -4 for each new parameter introduced into the fit (Bevington 1969)). The soft excess is 


2 This parameter is determined from Hi 21cm surveys. See http://heasarc.gsfc.nasa.gov/cgi- 
bin/Tools/w3nh/w3nh.pl for the column density calculator as well as Dickey & Lockman (1990). 
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typically parameterized by a blackbody component representing thermal emission in the central 
engine, but may also be modeled by bremsstrahlung emission or Comptonized emission from a 
thermal disk if either of these more complex forms gives a significantly better reduction in y} . The 
warm absorption is modeled using the same xSTAR-generated multiplicative table model described 
in Brenneman et al. (2007) and BR06 in the analyses of NGC 4593 and MCG-6-30-15. Solar 
abundances are assumed for all elements and the redshift is set at the source value. Some AGN 
show no need for a WA model, while others statistically require up to two physically separate WA 
components, each exhibiting a distinct column density (Ah) and ionization parameter (£,). 

Once the continuum has been properly modeled, we freeze the parameter values for all com- 
ponents except the power-law spectral index (F) and normalization (flux in phcm -2 s -1 ) and 
restrict our energy range from 2.5 — 10.0 keV in order to focus on the hard X-ray spectrum and the 
iron line region. Including the energies from 4.0 — 8.0 keV again, we check for residual emission 
(or absorption) lines from 6.4 — 6.97 keV, indicating the presence of unmodeled neutral and/or 
partly ionized Fe-Ka in the spectrum. Because these sources have been pre-determined to pos- 
sess significant, broadened iron emission, such residual emission features are seen in each case. 
We begin by attempting to fit the 6.4 keV line of neutral iron (and any other ionized lines) with 
a simple Gaussian feature, with the line width frozen at o = OkeV (i.e., intrinsically narrow) and 
the redshift again set to the source value. As expected, in each case we note significant residual 
“wings” remaining around the Fe-Ka feature after this narrow core was fit, indicating the presence 
of a broad component to the line. 


3.2. Fitting the Iron Line Profile 

To assess the robustness of the broad iron component and to obtain constraints on the BH 
spin, if possible, we follow our continuum fitting procedure (including modeling the narrow core 
of the iron line) by running each source through an “analysis tree” of progressively more complex 
modeling in the broad iron line region. Our procedure is outlined below: 

1. We begin by adding in a broad iron line described by the la or model (spin hardwired at 
a — 0.998) at the energy of the cold Fe-Ka line and refit, noting the corresponding change 
in x 2 /dof to verify whether a broad line is statistically warranted in the model. We also note 
the inner radius of disk emission: a broad line with the potential to diagnose BH spin should 
have an inner disk radius of r m j n < 6 r g . We relax this restriction here to allow for differences 
in fit that may be achieved through different models, requiring that r m j n < 20 r g . 

• If the broad line is not statistically warranted, we simply constrain the upper limit of 
the equivalent width of the line. 
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• If the broad line is at least marginally statistically warranted, we replace the la or 
component with a kerrdisk component (a = 0 — 0.998) and refit, noting again the 
corresponding change in % 2 /dof. The inner disk radius and BH spin as well as their 
statistical constraints should also be noted. 

2. Rather than simply modeling the Fe-Ka line in isolation, we replace the broad line compo- 
nent with a static ionized disk reflection spectrum (Ross & Fabian 2005) and fit. 

3. Starting from the static ionized disk reflection sprctrum, we add in the effects of relativistic 
smearing from the spacetime around a maximally rotating BH using the kdblur convolu- 
tion model ( a = 0.998). Again we note the statistical difference in the goodness-of-fit to 
determine the importance of relativistic effects (our criterion again being F-test statistics), as 
well as the inner radius of disk emission. This procedure gives us a rigorous measure of the 
degree to which broadening is required even when reflection is accounted for. 

• If the inner radius of disk emission is sufficiently close to the event horizon (again 
using r m i n < 20 r g ), the line may be sufficiently broad and relativistic effects may be 
sufficiently important to allow us to constrain BH spin. In this case, we replace the 
kdblur component with kerrconv (a = 0 — 0.998) in order to obtain statistical con- 
straints on rmin and a. 

We force the outer radius of disk emission to be r max = 400 r g (or 400 r ms in the case of 
kerrdisk and kerrconv, which parameterize disk emission radius in units of the radius of marginal 
stability) in all the broad line fits. The value itself is somewhat arbitrary; as long as a > 2 emission 
from this part of the disk is negligible. For the ionized disk component, we begin with the assump- 
tion of solar iron abundance and relative neutrality (Fe/solar = 1 and \ = 30 erg cm -1 s -1 ). We 
then relax these assumptions and allow these two parameters to fit freely, provided that reasonable 
statistical error bars can still be obtained for each parameter. For kerrdisk and kerrconv, we 
initially assume that the inner and outer portions of the disk emit under the same emissivity index: 
cxi = CC 2 , where the disk radiates as r~ a at any given radius. Again, if the fit is sufficiently robust, 
these indices are ’’untied” and allowed to each fit freely, but due to limitations in the number of 
photons we have for each data set, the only source with enough counts to support a scenario where 
<Xi ^ 0 C 2 is MCG— 6-30-15. 


4. Sample Results 


We have compiled the relevant parameter constraints from all of our spectral fits to nine Sy-1 
sources: MCG-6-30-15, MCG-5-23-16, Mrk 766, NGC 3783, NGC 4051, Ark 120, Fairall 9, 
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NGC 2992, and 3c273. The constraints on the relativistic disk parameters for each object are 
presented in Tables 2-4. Note that in each case a model including reflection from an ionized disk 
convolved with a relativistic smearing kernel provides the best statistical fit, and that in almost 
every case the kerrconv model (with arbitrary spin) provides a better fit than the kdblur model 
(spin fixed at a = 0.998). 

Previous studies have indicated that broad iron lines may be present in up to 42% of AGN 
observed with the XMM-Newton/EPIC-pn camera that have > 10 4 counts (Guainazzi et al. 2006). 
To assess whether evidence exists for a relativistically broadened iron line in each of our data 
sets, we first considered the laor fits in which BH spin is fixed at a = 0.998. While each source 
demonstrated a significant improvement in its global fit as compared to a model fitting only a 
narrow iron line core with a Gaussian component, the important parameter to evaluate in this case 
is the inner radius of disk emission. Roughly speaking, as stated in §3.2, if r m j n < 20 r g we can 
say with some confidence that there is substantial emission from the inner part of the accretion 
disk where relativistic effects (such as BH spin) become important in shaping the overall iron 
line profile. Out of our nine sources, four did not meet this criterion outright (MCG-5-23-16, 
NGC 3783, Ark 120, and Fairall 9). Allowing for arbitrary BH spin via kerrdisk instead of laor 
does help the situation somewhat: only NGC 3783 did not meet the r m j„ criterion in this fit, though 
for none of the sources is the global fit greatly improved over its laor value. Fairall 9 and 3c273 
have spins that cannot be constrained with the kerrdisk model. The laor and kerrdisk fits are 
presented in Table 2. 

We have also examined the question of how robust the presence of a broad line is in the 
spectrum when reflection is included. Beginning with a base continuum model including a static 
ionized disk reflection spectrum (Ross & Fabian 2005), we noted the overall goodness-of-fit and 
the residual features remaining, especially around the iron line region. We then convolved this 
model with relativistic effects, first using the kdblur smearing kernel (Table 3), then substituting 
in kerrconv (Table 4). A significant improvement in the global fit with kdblur indicates that rel- 
ativistic smearing is important in the system, and further significant improvement with kerrconv 
instead of kdblur indicates that although relativity is important, in the best fit scenario the BH spin 
deviates from its maximal value. Whenever possible we allowed the iron abundance and ionization 
parameter to fit freely in the ionized disk model, though in some cases these additional degrees of 
freedom prevented us from obtaining a reliable fit. For these sources it was necessary to assume a 
neutral disk (£, = 30 erg cm -1 s~' ) and solar iron abundance (Fe/solar=l): Ark 120, Fairall 9 and 
3c273. The kdblur and kerrconv fits returned similar results in terms of parameter constraints 
and overall fit to the iron line profile in each source, and in nearly all cases also produced a statisti- 
cally significant improvement in the overall goodness-of-fit, with the exception being NGC 3783, 
in which only a marginal improvement in the goodness-of-fit was seen. These results demonstrate 
that relativistic smearing is strongly required in nearly all of our sources, even after static reflection 
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is included. 


The ref lion model does not include a formal parameter for the reflection fraction of the 
disk, so in order to facilitate comparison between our results and those from studies using other 
models to parameterize reflection (e.g., pexrav), we have estimated i? re fl in our best-fit model 
(kerrconv( ref 1 )). To perform this calculation, we have assumed that for F re fl = 1, half of the 
radiation from the hard X-ray power-law is received and reflected by the disk, with the other half 
traveling directly to the observer. Reflection fractions of greater than unity then indicate that more 
than half of the power-law radiation illuminates the disk, implying anisotropic emission from the 
power-law source or perhaps strong light-bending of this emission near the BH. R rt n < 1 means 
that some of the power-law emission is either absorbed or scattered away before reaching the disk. 
We compute F re fl using the following formula: 


M'efl ^refl / ^ \ 

Apo Fpo \30 J 


( 1 ) 


Here /V re n and Npo denote the normalizations of the ref lion and power-law components from 
the best-fit model from 2.5 — 10.0 keV, (£/30) is the best-fit ionization parameter with respect to 
“neutral” in the ref lion component, and F re fl and F po are the fluxes contained in the ref lion and 
power-law components for unity normalization over the full wavelength range (0.001 — 1000 keV). 
The estimated reflection fraction for the best-fit model, F re fl, including errors, is shown in Table 4. 

We have obtained formal constraints on BH spin in four of the eight new cases: Ark 120 
(a = 0.63 — 0.68), Mrk 766 (a > 0.85), NGC 4051 ( a > 0.67) and 3c273 (a > 0.72). In comparison, 
the best fit to MCG-6-30-15 yields a > 0.96 from 2.5 — 10.0 keV and a > 0.987 based on our 
0.6 — 10.0 keV fit discussed in BR06. The two MCG-6 results differ slightly due to the different 
energy ranges modeled in BR06 vs. this paper. We were unable to constrain BH spin in MCG-5- 
23-16, NGC 3783, Fairall 9 and NGC 2992 in the kerrconv fits. 


There are several potential reasons why one might not be able to obtain a good spin constraint 
for a given source. Firstly, a lack of sufficient photons in the data set is a severe limitation for 
obvious reasons. This is certainly the case in 3c273, for example, and also plays a role in the 
spectra of NGC 2992 and Fairall 9. Generally speaking, the more photons one has to work with, 
the better the fit constraints will be. 


Secondly, if the iron line is not sufficiently broad, the relativistic parameters will be more 
difficult to constrain. A narrow iron line could originate from far out in the disk, away from the 
spacetime where relativistic effects become imprinted on the line profile. This scenario would 
leave us with no way to know the spin of the BH. Even iron lines from the inner regions of disks 
around non-spinning BHs show the spectral signatures of relativity, as shown in BR06, so if the 
line originates reasonably close to a non-spinning BH, we should be able to constrain its spin. 
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Therefore, we should not be surprised that the sources for which we are unable to get spin con- 
straints all share a common characteristic: constraints on r m j n that fall outside the inner portion of 
the disk, i.e., well outside of 6 r g . 

Our methodology uses the degree of relativistic smearing (and principally the broadening of 
the iron line) to constrain black hole spin on the assumption that the inner radius of the iron line 
emitting region cannot be smaller than the ISCO. Even allowing for some emission to originate 
within the ISCO will not significantly alter the iron line profile (Reynolds & Fabian 2008). How- 
ever, one can readily envisage scenarios whereby the iron line emission can truncate at radii larger 
than the ISCO, for example, due to a very highly ionized inner disk or a transition to a radiatively- 
inefficient accretion flow. This leads to an important (approximate) degeneracy in our model fits; 
since our spin constraints are primarily driven by modeling the maximum redshift of the iron line 
(rather than higher order details of line profile), the possibility that the line emitting region can 
truncate outside of the ISCO implies that our derived spin constraints will be in the form of lower 
limits. This is indeed what we find in three out of our four cases (excluding MCG-6-30-15). In 
the case of Ark 120, the accretion disk appears to be rather edge-on and (within the context of our 
spectral model) strong light bending effects around this highly inclined disk allow the degeneracy 
between a and r m i„ to be broken. However, for more face-on disks, exceptional signal-to-noise 
data (exceeding that possible from XMM-Newton ) is required to break this degeneracy via details 
of the line profile. 


AGN 

a 

a 

in 

r min (r g ) 

X 2 /dof 

MCG-6-30-15 

2.84-2.95 

5.93 - 6.37, 2.27 - 2.46, 10.28 - 13.99 

0.998 

0.97-0.98 

30-32 

30-31 

1.71 -2.11 
1.68-2.67 

940/1002(0.94) 

907/999(0.91) 

MCG-5-23-16 

1.98-4.30 

2.27-2.82 


21-37 

38-39 

18-44 

32-65 

1509/1490(1.01) 

1514/1489(1.03) 

Mrk 766 

1.78-2.05 

2.14-2.33 

0.998 

0.39-0.62 

85-86 

77-78 

1.24-1.82 
4.41 -4.67 

1542/1311(1.18) 

1539/1310(1.17) 

NGC 3783 

9.23 - 10.00 

0.998 

>0.89 

< 18 
17-18 

1.24-1.80 

315-388 

1741/1489(1.17) 

1732/1488(1.16) 

NGC 4051 

2.45-3.63 

4.24-10.00 

0.998 

<0.87 

30-39 

32-41 

3.96-6.98 

3.63-5.67 

1480/1397(1.06) 

1479/1396(1.06) 

Fairall 9 

1.40-2.70 

1.79-3.50 

0.998 

>50 

52-69 

1.24-25 

3.86-59 

672/692(0.97) 

671/691(0.97) 

Ark 120 

1.65-1.98 

0.998 

<0.69 

52-69 

67-69 

> 150.00 
3.94-6.43 

1639/1422(1.15) 
1620/1421 (1.14) 

NGC 2992 

1.31-3.15 

3.90-10.00 

0.998 

>0.97 

37- 55 

38- 49 

3.07-4.27 

1.67-6.22 

1399/1310(1.07) 

1399/1309(1.07) 

3c273 

2.83-6.37 

0.998 

<67 

<35 

1.24-3.18 

1.64-3.99 

228/267(0.854) 

229/266(0.86) 


Table 2: Comparison of the spectral fitting results for the 2.5 - 10.0 keV spectra using the laor and kerrdisk models. All 
AGN in our sample that are thought to possess broad iron lines are presented, as well as the MCG-6-30-15 data from BR06 
re-fitted with these models from 2.5 — lOkeV, in keeping with our current sample. For each source the top row represents the 
laor range in parameter value while the bottom row represents the kerrdisk range in parameter value. Ranges are given at 
90% confidence. Note that MCG-6-30-15 is able to be fit with a broken power-law emissivity index, which is listed under a as 
at, a 2 , ri>r (r g ). For all other sources, (Xi = (X 2 and rt, r = 6.0 r g . 











AGN 

Fe/solar 

^rcfl (ergcm -1 s -1 ) 

a 



X 2 /dof 

MCG-6-30-15 


249-273 




1296/1003(1.29) 



<39 

8.89-10.00 

39-40 

4.47-4.94 

887/1000(0.87) 

MCG-5-23-16 


<34 




1555/1492(1.04) 



<31 

1.05-4.40 



1485/1488 (1.00) 

Mrk 766 

1.00 

1742-3886 




1534/1313(1.17) 


1.00 

413-672 

1.71-2.04 

73-81 

1.24-2.96 

1495/1310(1.14) 

NGC 3783 

0.18-0.34 

<53 




1534/1490(1.03) 


0.14-0.23 

<46 

— 

<35 

> 102 

1525/1487(1.03) 

NGC 4051 

1.21-4.27 

<57 




1506/1398(1.08) 


1.29-10.00 

< 122 

3.08-4.29 

32-38 

1.24-7.36 

1476/1395(1.06) 

Fairall 9 

1.00 

30.00 






1.00 

30.00 

2.47-10.00 


1.24-8.25 

665/692(0.96) 

Ark 120 

1.00 

30.00 




1717/1425(1.20) 


1.00 

30.00 

1.59-1.83 

82-86 

1.24-2.02 

1576/1422(1.11) 

NGC 2992 

1.00 

30.00 




1414/1313(1.08) 


1.00 

30.00 

1.50-4.17 

17-42 

— 

1405/1310(1.07) 

3c273 

1.00 

30.00 




235/269(0.87) 


1.00 

30.00 

2.45 - 10.00 

<60 

1.24-11 

228/266(0.86) 


Table 3: Comparison of the spectral fitting results for the 2.5 - lO.OkeV spectra using the ionized reflection model reflion 
and kdblur(refl), which convolves the reflection spectrum with a smearing kernel from a maximally-spinning BH. For 
each column of data, the top row represents the reflion range in parameter values while the bottom row represents the 
kdblur ( ref 1 ) range in parameter values. For the reflion component, the value of the incident power-law spectral index is 
set equal to that of the continuum. Ranges are given at 90% confidence. 
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5. Discussion and Conclusions 

This paper explores the use of broad iron emission lines to constrain black hole spin, and 
makes the first steps toward estimating the spins in a small and relatively homogeneous sample of 
AGN. Following our previous work on MCG-6-30-15 (BR06), this study of an additional eight 
prominent Sy-1 AGN observed to have prominent iron lines has provided us with some intrigu- 
ing data on this subject and has also put us in an excellent position to apply our spectral fitting 
techniques to many other sources in the future. 


5.1. Implications of our Sample Results 

Though our total sample of nine Sy-1 AGN is far too small to allow us to draw any robust 
conclusions about the distribution of BH spins in this group of objects as a whole, this work has 
allowed us to design a viable method for evaluating the presence of a broad iron line in the data and 
using it to place constraints on the BH spin with our new models, kerrdisk and kerrconv. Based 
on our data analysis we have arrived at several important conclusions that we hope will inform 
subsequent BH spin surveys: 

1 . Guainazzi et al. (2006) found that the number of broad iron lines observed in AGN increases 
dramatically when the data have more than ~ 10 4 * counts. Similarly, we find that it is neces- 
sary to have a certain minimum number of photons in order to have a spectrum of the quality 
needed to extract spin information. Based on our results, this number appears to be on the 
order of ~ 10 6 photons. This makes sense, given that spectral line fitting of any kind is a 
statistical process and without enough data points one cannot expect to model a line with 
any statistical certainty. 

2. A broad iron line is needed in order to constrain spin. Simply put, a narrow line cannot 
represent emission from the inner part of the disk; if it were, it would show signatures of 
Doppler shifting and relativistic effects that combine to give good diagnostic disk lines such 
as Fe-Ka their characteristic shapes. If a line does not originate in the inner disk, then 
relativity is not important to its morphology, and it cannot be used as a spin diagnostic 
because spin is a relativistic effect. 

3. Five out of nine sources from our broad Fe-K line sample (including MCG-6-30- 1 5) showed 
a significant improvement in their global goodness-of-fit when relativistic smearing was 
added to a static ionized disk reflection spectrum. 


4. Across our sample of radio-quiet Sy-1 sources, we see a mix of spin constraints from mod- 

erate to very high. We do NOT see any non-spinning sources among those with robust spin 
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constraints, although it is possible that those sources without good spin constraints may be 
slowly-spinning. This shows that radio-loudness is not dependent on spin alone: it must 
be determined by at least one other parameter. Magnetic fields and/or thick disks are 
likely candidates: the Blandford-Znajek process is thought to be critical in jet formation, 
and most sources with radio jets are hypothesized to emanate from systems with truncated 
thin disks (perhaps ADAFs or other types of radiatively-inefficient accretion flows). 

5. Unfortunately, nature may be working against us in terms of measuring BH spin in radio-loud 
sources using the iron line method. Jet activity tends to overwhelm the spectrum along lines 
of sight to the broad line region, effectively burying the broad line and other disk signatures 
in many sources. Though some recent observations have noted the presence of a broad iron 
line in radio galaxies (e.g. 3c273, this work; 3c 120, Kataoka et al. 2007), for most radio-loud 
AGN it may be necessary to rely on other methods of diagnosing spin. 

Our five AGN with robust spin constraints all show a > 0.3, which seems in conflict with the 
theory presented by King, Pringle & Hofmann (2008) that BHs fed by accretion should tend to 
have spins in the range of a = 0. 1 — 0.3 if the accretion events occur in random directions over the 
lifetime of the AGN. Moreover, the rapid spin seen in spiral-hosted MGC-6-30-15, in particular, 
also refutes the hypothesis that rapid rotators should be found in giant elliptical galaxies (King et al. 
2008). Rather, our results tend to agree more with the Volonteri (2005) simulations predicting that 
BHs should grow mostly by accretion and evolve toward higher spins. While we are not yet able 
to comment on the mechanism by which our BH spins evolved, our results strongly argue against 
a scenario which leads to the majority of Seyfert AGN having slowly-spinning BHs. Given our 
limited sample size, however, more data are clearly needed in order to address these inconsistencies 
between theories and observations. 

The last bullet point will prove to be an interesting area of research in the coming years, espe- 
cially with newer, more sensitive, broader bandwidth instruments such as Suzaku and simultane- 
ous X-ray/radio studies of broad line radio galaxies. The question of how BH spin correlates with 
radio-loudness still remains unanswered. Sikora et al. have examined the radio-loud/radio-quiet 
dichotomy in AGN and argue that BH spin is a critical factor in determining the radio-loudness 
of an AGN. The authors conclude that both spiral-hosted and elliptical-hosted AGN show radio- 
loudness increasing with decreasing Eddington ratio and that the large host-morphology-related 
difference between the radio-loudness reachable by AGN in disk vs. elliptical galaxies can be ex- 
plained by the dual postulates that (1) the spin of a BH determines the jet outflow power, and (2) 
AGN BHs can reach large spins only in early type galaxies following major mergers (Sikora et al. 
2007). The latter seems to conflict with the hypothesis that BHs most likely reach near-maximal 
spin values through steady accretion since mergers can theoretically occur with random BH spin 
directions (Volonteri et al. 2005). Also, the extremely high spin observed in MCG-6-30-15, in par- 
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ticular, serves as evidence that high BH spins can be found in AGN hosted by disk-type galaxies. 
If our 3c273 result is found to be robust, we now also have an example of a radio-loud quasar with 
a moderately high BH spin. These conflicts between data and theory highlight the growing need 
for observations of BH spin in a large sample of AGN with varying physical characteristics (e.g., 
mass, radio-loudness, etc.). 


5.2. Formal Constraints on the Spin Distribution Function of Supermassive Black Holes 

Armed with our spin constraints, we are able to observationally assess the spin distribution 
function (SDF) of supermassive black holes in local AGN for the first time. In doing this, we must 
assume that the nine sources in our sample are an unbiased subset (in terms of black hole spin) 
of the parent population. Since they were selected on the basis of possessing known broad iron 
lines, this assumption amounts to saying that the detectability of broad iron lines is unaffected by 
black hole spin. Significantly larger samples are required in order to actually test this assumption. 
However, one selection effect that might be expected results from the decreased detectability of 
broad iron lines due to the blending of the iron line with the continuum; this effect would be 
expected to bias our sample towards more slowly spinning black holes (and hence narrower iron 
lines). 

Neglecting such selection effects, we proceed by parametrizing the SDF as a power (with 
index Q of the dimensionless spin parameter, 

/(«)=« + l)J- «>-l). (2) 

We employ a Monte Carlo technique to assess the probability distribution of C, given our data. For 
a given value of we conduct a trial by randomly drawing nine spins from the distribution in order 
to mimick the nine sources in our sample. The trial is deemed acceptable if the five spin constraints 
reported in Table 4 can be accommodated by the trial sample. This procedure is repeated 10 6 times 
for each value of £ in a grid spanning the range C, <E (—1,5) with spacing AC, = 0.01. Assuming 
a uniform prior for C, in this range, we then use the fraction of successful trials as a function of C, 
to construct the probability distribution for C, shown in Fig. 1 . This distribution has a broad peak 
around C, ~ 1 , with 90% confidence limits on C, of —0.01 < C, < 3.06. 

Being our only constraint that is not in the form of a lower limit, the spin measurement in 
Ark 120 is extremely important in setting the form of the C, probability distribution shown in Fig. 1 . 
Given this fact, it is interesting to examine constraints on the SDF if we exclude the Ark 120 mea- 
surement, e.g., allowing for (future) alternative explanations of the unusual Fe-K band structure in 
the X-ray spectrum of Ark 120. If we only have lower limits, the lack of any physical upper limit to 
the prior for C, requires us to phrase our statistical question differently. Instead of constructing the 
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OL 



0 2 4 


Fig. 1. — Probability distribution for the index C, of the spin-distribution function (f(a) = (£ + 
l)fl£). See text for details and caveats of how this probability distribution is contructed. 
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probability distribution for £, we instead use our Monte Carlo technique to examine the probability 
P\ 'r of obtaining a trial sample consistent with our observational results as a function of Given 
that (without Ark 120) our observational results are only in the form of lower limits, Pr — * 1 as 
C, — » °°. We find that P^ = 0.05 when £ = —0.30 and P £ = 0.10 when £ = —0.09. Thus, in this 
sense, our 90% and 95% confidence limits on C, are £ > —0.09 and C, > —0.3, respectively. 


5.3. Future Work 

Because we know so little about the spin distributions of black holes, much work remains to 
be done in this field. While the development of free-spin models such as kerrdisk and kerrconv 
has allowed us to begin attacking this question quantitatively and robustly for the first time among 
AGN sources, our survey represents only a first step in a much larger picture. Clearly, we need 
more sources with broad iron lines to form a larger statistical sample. We need better quality 
data with more photons and higher spectral resolution in order to better constrain the BH/disk 
parameters in broad iron line sources. We can get this through longer observations with active 
observatories such as XMM-Newton, Chandra and Suzaku, as well as coordinated observations in 
other wavebands (e.g., radio and optical) which can help us form a more complete picture of the 
physics at work in these systems. Suzaku, in particular, is already proving to be an invaluable 
tool in the X-ray astronomer’s arsenal with its ability to probe energies from 0.2 — 600 keV. This 
enables the examination of the soft spectrum, iron line region and reflection peak simultaneously, 
greatly aiding our ability to constrain the physical parameters of the accretion disk and isolate the 
effect of BH spin. 

We also need to expand our search radius to include stellar-mass BHs. These objects are 
thought to have their natal spins intact, so it will be interesting to see how the spin distribution in 
this population compares to that of AGN, especially given the physical similarities between AGN 
and GBHs in so many other areas. The advantage of the iron line method is its lack of dependence 
on BH mass, so this is an ideal technique to employ for such a comparison. Given the presence 
of other methods to diagnose BH spin, it will also be important to see how our results compare to 
those obtained using thermal continuum fitting, polarimetry and QPOs. Ideally, the spin constraints 
calculated using each different method should be consistent. 

Finally, we can improve the accuracy of the kerrdisk and kerrconv models by including 
the effects of relativity on emission produced from within the ISCO of the disk. In so doing, 
however, it is crucial to properly account for the ionization level and optical depth of the material 
in the plunging region. Because this gas is subject to greater tidal forces and incident blueshifted 
radiation than the gas outside the ISCO, it is likely to be quite rarefied, optically thin and highly 
ionized. Proper modeling of the gas physics in this region is necessary in order to evaluate the 
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overall contribution of emission from this region to the broad iron line. 
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A. Appendix: Source-by-source Model Fitting 

As mentioned in the previous Section, we have selected the AGN for our sample primarily 
based on the research performed by other groups on the robustness of broad iron line observations 
in various AGN and the likelihood of obtaining viable spin constraints from these sources using 
the iron line method (Nandra et al. 2007; Miller 2007). In particular. Miller (2007) discusses recent 
results from 30 Seyfert AGN in which relativistic disk lines have been detected and published, and 
divides those sources up into three “tiers” based on the nature and robustness of the detections. 
Many of the sources overlap with those presented in Nandra et al. Though we have examined 
observed spectra from all of the 9 AGN listed in Tier 1 (the most robust cases), NGC 3516 and 
NGC 4151 were not included in our study due to their extremely complex soft spectra (Turner 
et al. 2005; Schurch et al. 2004). We did analyze the other 7 sources in full: 3cl20, MCG-6-30- 
15, MCG-5-23-16, NGC 2992, NGC 4051, NGC 3783 and Mrk 766. 3cl20 is not included in 
this work, however, as our spectral fitting results for this object were inconclusive in several ways 
due to in part to a paucity of photons. We have also examined 3 other sources from Tier 3 that are 
mentioned prominently in Nandra et al. : Fairall 9, Ark 120 and 3c273. 

We present the results of our spectral analysis for each source in its own subsection below. 
Note that the reflection fractions quoted are not fitted values, but rather are estimates calculated 
using the method described in §4 for the purpose of comparison to other studies which may model 
disk reflection differently (e.g., with pexrav). Detailed spectral fitting results for MCG-6-30-15 
can be found in BR06 and are therefore not included herein. 

A table of best-fit parameters for each source is included, as are figures representing the 
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residual broad Fe-K line feature, the best fit to the spectrum, and the best-fit model components 
for each object. The line residuals represent the data-to-model ratio for the best-fit model after 
removal of the kerrconv( ref 1 ) components used to fit the broad Fe-K line; the narrow line core 
is left in the model. 


A.0.1. MCG-5-23-16 

MCG-5-23-16 is a moderately absorbed Seyfert galaxy of intermediate type (Sy-1.9). It is 
relatively nearby at a redshift of z = 0.0085, and with a typical 2 — lOkeV flux of F 2-10 ~ 8 x 
10 -11 ergcm _2 s _1 it is one of the brightest known Seyfert galaxies (Reeves et al. 2006). The 
source has been observed previously to have a soft excess below ~ 1 keV and an absorbing column 
of Ah ~ 10 22 cm -2 (Dewangan et al. 2003; Balestra et al. 2004), and ASCA observations have 
indicated the presence of a broad Fe-Ka line of EW ~ 200 eV (Weaver et al. 1997, 1998). The 
line was successfully modeled with a narrow core at the rest-frame energy of 6.4 keV and a broad 
component superposed on it. This feature was modeled with similar success in the XMM-Newton 
observations of Dewangen et al. and Reeves et al. , as well as the Suzaku observation also published 
by these authors (Reeves et al. 2007). 

In December 2005, MCG-5-23-16 was simultaneously observed with XMM-Newton, Chan- 
dra, Suzaku and RXTE. The XMM-Newton results are reported by Reeves et al. , and our data 
reduction followed that work (Reeves et al. 2006). The EPIC-pn instrument had a net exposure of 
96 ks and returned ~ 2.2 x 10 6 photons after filtering. 

As suggested by Reeves et al. , a simple photoabsorbed power-law fit does not model the 
continuum well, especially below ~ 1 keV. Adding in a source of soft thermal emission also leaves 
prominent residuals and does not adequately account for the shape of the soft excess, so following 
the lead of Reeves et al. we have employed a two-power-law model (Reeves et al. 2006): one 
of the power-law components is affected only by Galactic photoabsorption (Ah = 8 x 10 20 cm -2 
as per Reeves et al. ), and hence leaves the AGN system effectively unscattered, while the other 
component is also subject to absorption intrinsic to the AGN, in this case with a fitted absorbing 
column density of Ah = 1.19 x 10 22 cm -2 . This component also experiences scattering within 
the system. We set both power-law photon indices to be equal, indicating that the two power-law 
components originate from the same physical reservoir of photons (with T = 1 .66). Calculating the 
ratio of the flux of the scattered power-law component to the unscattered power-law component 
(from 0.6 — 10.0 keV, fitting only the continuum) yields an estimate of the optical depth of the 
scattering plasma. In this case x = 7.37 x 10~ 4 phcm -2 s~ 1 / 1 .83 x 10~ 2 phcm~ 2 s _1 « 0.04, so 
the scattered fraction is low, implying that the electron plasma is optically thin. 
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The spectrum also shows evidence for two prominent, narrow emission lines of iron Ka and 
Kp (at E Ka = 6.4 keV and E K?1 = 7.0 keV, with EW Ka = 167 eV and EW Kf> = 1 04 eV, respec- 
tively), as well as three narrow absorption lines of intermediately ionized iron at E\ = 7.24 keV, 
E 2 = 7.48 keV and £ 3 = 7.85keV {EW X = -79eV, EW 2 = -121 eV, and EW 3 = -143eV, re- 
spectively). 

After successfully modeling the continuum and narrow line parameters for MCG-5-23-16, 
we restricted our attention to the hard spectrum (2.5 — 10.0 keV) and analyzed the broad line com- 
ponent of this source using the method outlined above in §3.2. Note that eliminating soft energies 
from the fit removed the statistical need for a second, soft power-law in the fit. For a full listing 
of best-fit model parameters and error bars for MCG-5-23-16, see Table 5. The residual iron line 
feature, best-fit model and best-fit model components for the hard spectrum are shown in Fig. 2. 
The best-fitting model for MCG-5-23-16 is the ionized disk reflection spectrum (Ross & Fabian 
2005) convolved with kerrconv. Although we were not able to constrain the BH spin in this 
source, we did achieve a fit constraint on the inner radius of the disk emission: r m i n ~ 16r ms . This 
corresponds to a range of ~ 36 — 96 r g for a BH spin between a = 0.0 — 0.998. Given that the 
effective inner edge of the disk is not within the radius of marginal stability for a Schwarzschild 
BH, it is not surprising that we were not able to constrain the value of the BH spin in this source. 


A.0.2. NGC3783 

NGC 3783 is a bright, nearby Sy-1 galaxy at a redshift of z = 0.0097. It was first detected 
in X-rays with the Ariel-V all-sky survey (McHardy et al. 1981) and subsequently in the high 
Galactic latitude survey conducted by HEAO-1 (Piccinotti et al. 1982). Since these early detec- 
tions, there have been many observations of NGC 3783 with higher resolution instruments: A 
ROSAT observation of the source showed evidence of an ionized absorber in the soft band (Turner 
et al. 1993), which was confirmed during ASCA observations (George et al. 1995, 1998). High- 
resolution grating observations of NGC 3783 with Chandra and XMM-Newton followed, unveiling 
the soft spectrum of this source in detail (Kaspi et al. 2000, 2001, 2002; Blustin et al. 2002; Behar 
et al. 2003). The higher signal-to-noise of XMM-Newton, in particular, has also enabled the iron 
line to be studied extensively in this source (Reeves et al. 2004). Using two observations from 
December 2001, Reeves et al. have obtained ~ 240 ks of data on NGC 3783. Their global fits to 
the merged EPIC-pn spectrum are currently second in length and depth only to the ~ 330 ks ob- 
servation of MCG-6-30-15 (Fabian et al. 2002). In these observations, the source has an average 
2 — lOkeV flux of F 2-.10 = 6.8 x 10 -11 erg cm -2 s -1 . The spectrum is noted to have iron line 
peaks at 6.4 keV and 7.0 keV, representing neutral Fe-Ka (EW ~ 120eV) and a blend of ionized 
Fe-Ka and Fe-K(3 {EW ~ 35 eV), respectively. A strong absorption line of highly ionized iron is 
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Fig. 2. — Top left: The 2.5 — lO.OkeV spectrum of MCG-5-23-16 fit with a two-power-law model 
modified by Galactic photoabsorption; data-to-model ratio. Top right: Best-fit model for MCG-5- 
23-16, including our continuum model, two narrow emission lines, and five narrow absorption lines 
as well as an ionized disk reflection spectrum convolved with our kerrconv relativistic smearing 
kernel. Bottom: A vF v plot depicting the relative flux in each of the best-fit model components in 
the spectrum. 
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Table 5: Best-fitting model parameters for the 2.5 - lO.OkeV spectrum of MCG-5-23-16, including 
components and parameter values from 0.6 - 1.5 keV for completeness. The energies from 1.5 - 
2.5 keV were not included due to the presence of absorption edges from the XMM-Newton mirrors. 
Error bars are quoted at 90% confidence. For the z gauss lines, we required each to be intrinsically 
narrow, i.e., o = 0.0. Redshifts were frozen at the cosmological value for the source, in this case 
z = 0.0085. Note the absence of warm absorption and a soft excess in this object. When no error 
bars are quoted, the parameter in question is frozen at the given value. 
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also observed at 6.67 keV that exhibits variability in direct correlation with that of the continuum 
flux. A weak red wing of the 6.4 keV line is noted by the authors, but once warm absorption is 
taken into account in the system the requirements for a broad Fe-Ka component are significantly 
reduced (Reeves et al. 2004). 

We analyzed the first of the two December 2001 XMM-Newton observations. We chose not 
to merge the event files of both data sets to avoid the uncertainties inherent in so doing, instead 
focusing on only one data set with an effective exposure time of ~ 1 20 ks and a total of ~ 1.8 x 10 6 
photons after filtering. This number of counts is still large enough to obtain valid statistical fits to 
the model parameters used. We follow the data reduction steps of Reeves et al. , excluding the last 
part of the observation due to contamination from a background flare. Our analysis follows in the 
same spirit as these authors, but in order to remain consistent with the method described above in 
§3 we do not make any assumptions about the soft spectrum and continuum emission based on 
prior RGS results (Blustin et al. 2002). 

We begin our analysis of the time-averaged spectrum with a simple photoabsorbed power-law 
fit using the Galactic absorbing column of Ah = 8.5 x 10 20 cm -2 , and neglecting the energy ranges 
associated with the mirror edges and the iron lines, as in §3. This fit leaves obvious residuals in 
all parts of the spectrum, however, and is especially poor for the soft energies below ~ 1.5 keV. 
Continuum curvature associated with a soft excess and warm absorption are both clearly evident. 
As discussed by Reeves et al. , the fit is much improved with the addition of a thermal blackbody 
component (kT ~ 0.07 keV, as compared with their value: kT ~ 0.09 keV). We approach the 
question of the warm absorption similarly as well, using our XSTAR table model to parameterize 
the column density and ionization level of the absorbing medium, as do Reeves et al. As with 
these authors, we find statistical evidence for a two-zone warm absorption structure, though our 
values for the column densities and ionization parameters of these zones vary significantly from 
those of Reeves et al. We find N\\\ = 6.84 x 10 23 and log^i < 0.01, with Ah2 = 1-58 x 10 23 
and log ^2 = 1-40, whereas Reeves et al. find Ahi = 6.00 x 10 20 and log^i = 0.30, with Ah2 = 
4.60 x 10 22 and log ^2 = 2.90. These differences are not surprising, however, because Reeves et 
al. base their fits on the RGS results for this source (Blustin et al. 2002). We do not use any a 
priori information to augment or guide the EPIC-pn spectral fits due to calibration uncertainties. 
Finally, we allowed the value of the neutral hydrogen absorbing column to vary in an effort to 
improve the goodness-of-fit and found that the spectrum preferred a higher value than the Galactic 
column: Ah = 3.60 x 10 21 cm -2 . This was also in contrast to the Reeves et al. result, where the 
cold hydrogen column was held at the Galactic value and not permitted to vary. 

Reeves et al. found two emission peaks in the spectrum at 6.4 keV (EW ~ 120 eV) and 7.0 keV 
(EW ~ 35 eV), representing cold Fe-Ka and likely a blend of ionized Fe-Ka and Fe-K(5, respec- 
tively. The latter line, in particular, showed no appreciable variance over the observation, leading 
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the authors to postulate a relatively distant origin for the lines away from the central parts of the 
accretion disk. A 6.67 keV absorption feature of highly ionized iron was also seen which did 
appear to vary with time and was strongest when the continuum flux was highest, suggesting an 
origin in the region of the warm absorber within 0. 1 parsec of the nucleus. Though inclusion of 
the WA and this feature did lessen the statistical case for a broad iron line in NGC 3783, Reeves 
et al. nonetheless identified a residual broad feature which they fit using a diskline model with 
a ~ 3.3, i ~ 19°, r m j n = 6.0 r g , and EW ~ 58 eV. 

We also detected the emission and absorption features discussed by Reeves et al. , though 
we found that the equivalent widths for the Fe-Ka, Fe-Ka/KP blend and ionized iron absorption 
lines differed from the fits performed by these authors: EWko. = 79.0 eV, EW Ka / K $ = 23.9 eV and 
£W a bs = — 25.1 eV. We did not test the time variance of the absorption component because such 
an examination is beyond the scope of this work. It should also be noted that, contrary to the 
Reeves et al. analysis, we required the core of the 6.4 keV line to be intrinsically narrow, so it is not 
unusual that we obtain a smaller equivalent width than Reeves et al. When we included the broad 
component of the neutral Ka line, the fit improved dramatically. Our best fit to the hard spectrum 
of NGC 3783 convolved an ionized disk reflection spectrum with kerrconv, as can be seen in 
Table 6 below, but due to the relative narrowness of the iron line we were not able to constrain the 
BH spin in this source. Not surprisingly, our constraints on the radial extent of the disk showed that 
/"min < 390 r ms , or < 985 r g . A line emitted from so far out in the disk would not exhibit relativistic 
signatures. 


A.0.3. Mrk 766 

Mrk 766 is designated as a classic, bright NLS 1 galaxy (Sy-1 .5) with a redshift of z = 0.0129 
and a typical 2 — lOkeV flux of Z^-io ~ 2.5 x 10 -11 erg cm -2 s -1 (Pounds et al. 2003). Previous 
X-ray observations of this source have provided contradictory evidence on the detection of a broad 
Fe-Ka feature. Although not one of the most convincing cases, Mrk 766 was included in an ASCA 
spectral survey of bright Seyfert galaxies showing evidence of relativistic iron lines (Nandra et al. 
1997). A separate analysis of simultaneous ROSAT and ASCA observations (Leighly et al. 1996), 
however, showed the X-ray spectrum to be described by a power-law of index increasing strongly 
with flux from T ~ 1.6, but with only a narrow Fe-Ka emission line (EW ~ 100 eV at 6.4 keV). A 
later observation with BeppoSAX found a steeper power-law (T ~ 2.2), and evidence for an absorp- 
tion edge at ~ 7.4 keV (Matt et al. 2000), implying strong reflection from intermediately ionized 
material. Interestingly, based on XMM-Nexvton observations of the source. Turner et al. found that 
energy-time maps of Mrk 766 reveal a periodic energy shift in an ionized component of Fe-Ka 
emission, with a period of ~ 165 ks. This can be interpreted as evidence for emission from orbit- 
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Fig. 3. — Top left: The 2.5 — 10.0 keV spectrum of NGC 3783 fit with a power-law model modified 
by Galactic photoabsorption, along with a two-zone warm absorber and a blackbody soft excess. 
Top right : The 2.5 — lO.OkeV best-fit model for NGC 3783, including our continuum model, two 
narrow emission lines and a narrow absorption line as well as an ionized disk reflection spectrum 
convolved with our kerreonv relativistic smearing kernel. Bottom: A vF v plot depicting the 
relative flux in each of the best-fit model components in the spectrum. 
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Table 6: Best-fitting model parameters for the 2.5 - lO.OkeV spectrum of NGC 3783, including 
components and parameter values from 0.6 - 1.5 keV for completeness. The energies from 1.5 - 
2.5 keV were not included due to the presence of absorption edges from the XMM-Newton mirrors. 
Error bars are quoted at 90% confidence. For the zgauss lines, we required each to be intrinsically 
narrow, i.e., a = 0.0. Redshifts were frozen at the cosmological value for the source, in this case 
z = 0.0097. Note the two-zone warm absorber and blackbody soft excess present in this object, 
similar to many other Sy-1 sources. 
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ing gas within ~ 100 r g of the central BH. A likely explanation is that this gas represents a hot spot 
on the disk illuminated by magnetic reconnection (Turner et al. 2006). 

We examine the May 2001 XMM-Newton observation of Mrk 766 taken by Mason et al. , 
from which the RGS results were published in 2003 (Mason et al. 2003) and the EPIC spectra 
were explored more thoroughly by other authors (Pounds et al. 2003; Turner et al. 2006). In our 
re-analysis, the EPIC-pn data have an effective exposure time of 128 ks and the filtered data yield 
~ 2.92 x 10 6 photons. 

Mrk 766 bears many spectral similarities to MCG-6-30-15, showing significant complexity 
beyond a simple photoabsorbed power-law fit. There is clear evidence for a thermal (bbody: kT = 
0.08 keV) soft excess below ~ 1 keV as well as two distinct physical zones of low (Ahi = 5.26 x 
10 22 cm -2 , log^i = 0.21) to moderately ionized (Ah2 = 9.50 x 10 23 cm -2 , log £2 = 3.04) intrinsic 
absorption. These components exist in addition to cold absorption by the Milky Way (Ah = 1 .7 1 x 
10 20 cm -2 ). The power-law in this source has a spectral index of T = 2.75, comparable with the 
result published by Pounds et al. Unfortunately, these authors did not probe the spectrum below 
3 keV, so we cannot compare our soft excess or warm absorption parameters to theirs. 

The narrow Fe-Ka core of Mrk 766 has a rest-frame energy of 6.4 keV, as expected, and 
an EW = 33.7 eV, again comparable with the findings of Pounds et al. Fitting this line with a 
Gaussian component leaves significant residuals strongly indicative of the presence of a broad 
line. This broad component is most successfully fit with a kerrconv(refl) model, yielding 
comparable results to the disk reflection model employed to model the broad iron line in Pounds 
et al. , though with different fit parameters. Our model also detects no absorption features above 
8 keV, as discussed by Pounds et al. , possibly due to advances in modeling the reflection spectrum. 
The BH spin, to 90% confidence, is relatively high at a > 0.85 and the inner edge of the accretion 
disk is constrained to r m j n < 2.25 r ms . The best-fit results are shown in Table 7. Fig. 4 shows 
the iron line residual, best fit to this residual, and relative contributions of the individual model 
components, respectively. 


A.0.4. 3c273 

3c273 is the most distant source in our sample at a redshift of z — 0.1583. Classified as 
a bright variable quasar, this AGN displays a strong jet and during epochs of radio-loudness it 
possesses the flat X-ray spectrum of a blazar, with highly beamed jet emission (Tiirler et al. 2006). 
When the jet is reduced in strength, however, this source has been observed to exhibit Sy-l-type 
accretion disk signatures such as a broad iron line. A minimum state occurred in March 1986 and 
allowed Robson et al. (1986) to identify a new near-infrared spectral component. An even better 
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Fig. 4. — Top left: The 2.5 — 10.0 keV spectrum of Mrk 766 fit with a continuum composed of a 
power-law model modified by Galactic photoabsorption, a two-zone warm absorber model and a 
blackbody soft excess. Top right: The 2.5 — 10.0 keV best-fit model for Mrk 766, including our 
continuum model and an ionized disk reflection spectrum convolved with our kerrconv relativistic 
smearing kernel. Bottom: A vF v plot depicting the relative flux in each of the model components 
in the spectrum. 
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Table 7: Best-fitting model parameters for the 2.5 — 10.0 keV spectrum of Mrk 766, including 
components and parameter values from 0.6 — 1.5keV for completeness. The energies from 1.5 — 
2.5 keV were not included due to the presence of absorption edges from the XMM-Newton mirrors. 
Error bars are quoted at 90% confidence. We required the zgauss line representing the core of the 
Fe-Ka feature to be intrinsically narrow, i.e., c = 0.0. Redshifts were frozen at the cosmological 
value for the source, in this case z = 0.0130. A two- zone warm absorber and blackbody soft excess 
are present in this object, as in MCG-6-30-15 and many other Sy-1 sources. 
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opportunity arose in early 2004, when the sub-millimeter flux of 3c273 was observed to be almost 
two times lower than in 1986. This new minimum triggered a slew of simultaneous observations 
with instruments in all wavebands such as INTEGRAL, XMM-Newton and RXTE, among several 
other optical, radio and sub-millimeter telescopes. The 2— lOkeV flux during this period was 
F 2 _io = 6.7 x 10 -11 erg cm -2 s -1 . Such a low flux in 3c273 has only been measured twice in 
the past, by Ginga in July 1987 (Turner et al. 1990) and by BeppoSAX on 18 July 1996 (Haardt 
et al. 1998), the latter of which was coincident with the low sub-mm flux mentioned above. This 
X-ray/sub-mm correlation strongly supports a synchrotron self-Compton origin for the X-ray jet 
emission in 3c273. 

Tiirler et al. obtained a 20 ks XMM-Newton observation of 3c273 during this jet-minimum 
state in June 2004. We use their thin-filter EPIC-pn observation in the interest of collecting as 
many photons as possible, though due to photon pile-up it is necessary to exclude the centralmost 
region of the source, as detailed by the authors (Tiirler et al. 2006). In our filtered data set we 
capture ~ 7.5 x 10 5 photons. 

Tiirler et al. note the inadequacy of a simple photoabsorbed power-law fit to the data. Their 
best fit to the continuum is achieved using two power-law components: Thard = 1.63 ±0.02 and 
r s 0 f t = 2.69 ±0.06, with the hard component flux ~ 2.3 times the soft flux. Both components were 
modified by Galactic photoabsorption with N\\ = 1.79 x 10 20 cm -2 . We based our initial contin- 
uum fit on theirs, with two power-law components: Thard = 1. 721^28 anc * r s0 f t = 3.01 +q 2 j, both 
consistent with the Tiirler et al. results. The ratio we calculate between the fluxes of the hard and 
soft power-law components is 4.18, however. This is significantly higher than that calculated by 
Tiirler et al. , and the error bars on the spectral indices, in particular, led us to consider eliminating 
the soft component. Our best-fitting continuum model, therefore, has only one power-law compo- 
nent at T = 2.03 and exhibits the same statistical goodness-of-fit as the two power-law model. 

Evidence for excess emission from 2.5 — 7.0 keV was also found by Tiirler et al. , supporting 
the presence of a broad Fe-Ka line when 3c273 is in a jet-minimum state. The authors quote 
this excess as significant at the 6a level with an integrated flux of 2.6 ±0.4 x 10 -4 ph cm -2 s -1 , 
corresponding to an EW = 166±26eV. These values are consistent with those reported from 
previous observations of the source in a jet-minimum state (Kataoka et al. 2002; Page et al. 2004; 
Yaqoob & Serlemitsos 2000). The breadth of the excess is not satisfactorily fitted by a Gaussian 
line or a diskline component from a non-rotating BH. If the line extent is real, the only remaining 
explanation, according to the authors, is that it is emitted from around a near-extreme Kerr BH. 
The sharp edge of the iron line at 7 keV suggests that the angle of inclination of the accretion disk 
is 35 - 40° (Tiirler et al. 2006). 

We find evidence for a similar excess around the Fe-Ka line in our analysis. Interestingly, 
we do not see a narrow emission line at 6.4 keV, but rather the narrow line core seems to coincide 
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with the He-like line of Fe-Ka at 6.66 keV with an EW = 73.3 eV. This may mean that the gas 
in this system is highly ionized, which would not be surprising in such an active source. Visually, 
the residuals left over after fitting this narrow line suggest the presence of a broad component. 
Modeling this component with our analysis method provides a modest improvement in the global 
goodness-of-fit compared with other sources, perhaps due to the distance of the source and the 
paucity of the data. With a 20 ks observation we have collected < 3 x 10 5 photons from 2.5 — 
lOkeV. Given more observing time and a larger number of (pile-up free) photons, our line analysis 
would be significantly improved. Taking into account these caveats about the robustness of the 
broad line, our kerrconv(refl) model still provides the best fit, as shown in Table 8. Fig. 5 
shows the iron line residual, best fit to this residual, and relative contributions of the individual 
model components in the best fit, respectively. To 90% confidence, a > 0.72 and r m j n < 2.30 r ms , 
but longer, more sensitive observations are necessary in order to further substantiate these results. 


A.0.5. NGC2992 

NGC 2992 is a Sy-1.9 galaxy that appears to show a broad iron line even though it is highly 
obscured. This source has been the subject of intense study due to the variability of its X-ray 
emission (Gilli et al. 2000). In 1997 and 1998, BeppoSAX caught NGC 2992 transitioning from 
a Compton-thick to a Compton-thin state, resulting in an order of magnitude increase in its X-ray 
luminosity as well as a qualitative difference in its spectral appearance, with more disk features 
(e.g., broad lines) being seen in several wavelengths. The source has a redshift of z = 0.0077. In 
the two BeppoSAX pointings, NGC 2992 displays a 2 — lOkeV X-ray fluxes of 0.63 and 7.4 x 
10 -11 erg cm -2 s -1 (Gilli et al. 2000). 

We use the May 2003 XMM-Newton observation of NGC 2992, totaling 29 ks, which trans- 
lates to ~ 6.02 x 10 5 photons in the filtered data set. This observation appears to be unpublished, so 
all the fit values referenced herein are our own, compiled using the reduction and analysis methods 
detailed in §3. 

Upon first inspection, it is immediately clear that NGC 2992 is heavily absorbed below 
~ 2keV. Though statistical evidence exists for a soft excess, none of the telltale signatures 
of intrinsic warm absorption are present, and the flux decreases precipitously at soft energies. 
When a simple photoabsorbed power-law model plus blackbody emission is applied to the con- 
tinuum, the remaining residuals on the soft end indicate that some absorption remains unac- 
counted for, but allowing the column density of hydrogen to fit freely above the Galactic value 
of Nh — 5.26 x 10 20 cm -2 neatly corrects the discrepancy. The final continuum parameters are 
Nil = 5.11 x 10 21 cm -2 , r= 1.69 and kT = 0.05 keV with the power-law flux 1.4 times greater 
than that of the blackbody component. 



-35- 





Fig. 5. — Top left: The 2.5 — 10.0 keV spectrum of 3c273 fit with a power-law model modified by 
Galactic photoabsorption. Top right: The 2.5 — 10.0 keV best-fit model for 3c273, including our 
continuum model and an ionized disk reflection spectrum convolved with our kerrconv relativistic 
smearing kernel. Bottom: A vF v plot depicting the relative flux in each of the model components 
in the spectrum. 
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Model Component 

Parameter 

Value 

phabs 

N H (cm- 2 ) 

1.79ii»»xI0- 2 

po 

Tpo 

flux(phcm _2 s _1 ) 

3.82«;g x 10- 2 

zgauss 

E(keV) 

flux(phcm -2 s -1 ) 
EW(eV) 

6.66 

9.30+2 x 10- 5 

73.30+5^ 

kerrconv 

«i 

«2 

n>r(>ms) 

a 

in 

r min (^ms) 
r max ms) 

4 22+ 
^• zz -1.48 

4 22+ 1-78 
*- LL -\ .48 

6.0 

> 0.72 
<63 
<2.30 
400 

ref 1 

Fe/solar 

£refl (ergcm -1 s -1 ) 

Uefl 

^refl 

flux(phcm -2 s _1 ) 

1.0 

< 355.88 

2 03+ 0 30 

4 76 +28 - 34 
^• /0 -3.84 

U0t\% 2 x 10- 4 

% 2 /dof 


288/265(0.86) 


Table 8: Best-fitting model parameters for the 2.5 — 10.0 keV spectrum of 3c273, including compo- 
nents and parameter values from 0.6 —1.5 keV for completeness. The energies from 1.5 — 2.5 keV 
were not included due to the presence of absorption edges from the XMM-Newton mirrors. Error 
bars are quoted at 90% confidence. We required the zgauss line representing the core of the Fe- 
Ka feature to be intrinsically narrow, i.e., o = 0.0. Interestingly, this line was found at a moderate 
ionization state of iron (6.66 keV) rather than at the neutral rest-frame energy of 6.4 keV. Redshifts 
were frozen at the cosmological value for the source, in this case z — 0.1583. No warm absorption 
or soft excess is seen in this object. Note the large error bars on the value of R re f\ that we calculate. 
These are likely an artifact of the relatively low photon count for this observation. 


-37- 


The Fe-Ka line is present with a narrow core at 6.4 keV and EW — 38.8 eV. Significant 
residuals remain surrounding the core, however, and the best statistical fit is achieved with a 
kerrconv(refl) model. See Table 9 for model parameters and error bars. Fig. 6 shows the 
iron line residual, best fit to this residual, and relative contributions of the individual model com- 
ponents for the best fit. Note that the BH spin cannot be constrained. Constraints are also dif- 
ficult to achieve on the disk emissivity index a, especially, though we have constrained the ef- 
fective inner disk radius to r m ; n < 11.78r ms . These large error bars are seen due to the paucity 
of photons in the hard band (2.5 — lOkeV) in this observation: with only 29 ks and a flux of 
^ 2 . 5 -io = 8-06 x 10" 11 erg cm" 2 s" 1 , we have < 2 x 10 5 photons to use in spectral fitting. This is 
only a small fraction (< 20%) of the number of counts we have at our disposal for the same energy 
range in MCG-6-30-15, so it is not surprising that our statistics are not as sound in this case. 


A.0.6. NGC4051 

NGC 4051, like NGC 2992, is a heavily absorbed Seyfert AGN (Sy-1.5; NLS1) seen to vary 
significantly in flux over the course of several observations (Guainazzi et al. 1996; Lamer et al. 
2003). The source is at a redshift of z = 0.0023 and has a typical flux on the order of a few times 
10" 11 erg cm" 2 s" 1 , though as stated above, this flux can vary significantly on a variety of time 
scales, along with the spectral characteristics of the source. Unusually low flux states in this object 
can last for weeks to months, during which time the X-ray spectrum shows a hard continuum 
power-law of spectral slope T ~ 1, but is dominated by a softer component at lower energies with 
T ~ 3 (Uttley et al. 2004). A highly broadened and redshifted iron line has also been noted in the 
low flux state of NGC 4051 with RXTE, suggesting that reflection features from the accretion disk 
close to the BH may remain constant in this source in spite of the large variations in continuum 
properties (Uttley et al. 2003). 

NGC 4051 was observed with XMM-Newton in May 2001 for a duration of 1 17 ks, yielding 
~ 3.27 x 10 6 photons in the filtered data set. The EPIC-pn results were first reported by Mason et 
al. , and suggested a continuum described by a power-law “pivoting” around 100 keV, according to 
a simultaneous observation with RXTE. Ultraviolet emission from the Optical Monitor on XMM- 
Newton was found to lag the X-ray emission by ~ 0.2 days, indicating that it is likely reprocessed 
X-ray emission. The X-ray emission itself showed variability on time scales as small as 1 — 2 
hours (Mason et al. 2002). These results were expanded upon by Pounds et al. , who noted the 
intermediate flux of the source at this time and validated the presence of both an iron line and a 
thermal soft excess during the May 2001 observation (Pounds et al. 2004). Ponti et al. recently 
re-analyzed this observation as well as a lower-flux pointing from November 2002 and reinforced 
the veracity of this model, while also considering the comparable efficacy of a model dominated 
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Energy (keV) 


Fig. 6. — Top left: The 2.5 — 10.0 keV spectrum of NGC 2992 fit with a continuum composed 
of a power-law model modified by Galactic photoabsorption as well as a blackbody soft excess. 
Top right : The 2.5 — 10.0 keV best-fit model for NGC 2992, including our continuum model and 
an ionized disk reflection spectrum convolved with our kerrconv relativistic smearing kernel. 
Bottom: A vF v plot depicting the relative flux in each of the model components in the spectrum. 
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flux (ph cm -2 s -1 ) 
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30.0 

1 6 Q + 0-02 

l-o^_0.03 

0 64 +0 ‘ 22 
0.20 
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X 2 /dof 


1402/1312(1.07) 


Table 9: Best-fitting model parameters for the 2.5 — 10.0 keV spectrum of NGC 2992, includ- 
ing components and parameter values from 0.6— 1.5 keV for completeness. The energies from 
1.5 — 2.5 keV were not included due to the presence of absorption edges from the XMM-Newton 
mirrors. Error bars are quoted at 90% confidence. We required the 6.4 keV zgauss line core to be 
intrinsically narrow, i.e., a = 0.0. Redshifts were frozen at the cosmological value for the source, 
in this case z = 0.0077. No evidence for warm absorption is detected, though a soft excess is found 
in this object, as in many other Sy-1 sources. 
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by ionized reflection from radii quite close to the BH (Ponti et al. 2006). 

We have also focused on the May 200 1 observation of NGC 405 1 , following the reduction and 
analysis of Mason et al. and Ponti et al. but using modem calibration files. A two power-law model 
parameterized the continuum with much greater accuracy than either a single power-law, a broken 
power-law or a “pivoting” power-law of the type employed in previous X-ray spectral analyses. 
Hard and soft components were detected: Thard = 2.01 and r so f t = 5.20. The hard/soft flux ratio 
for these components is 1.21, and each is modified by photoabsorption from neutral hydrogen 
with A/fi = 7.63 x 10 20 cm~ 2 . This column density is greater than the Galactic value by a factor 
of ~ 6, indicating the presence of cold absorption within the NGC 4051 system. Additionally, 
there is statistical evidence to support the presence of a warm absorber within the AGN as well: 
a significant improvement in the global goodness-of-fit is seen with the inclusion of one of our 
XSTAR multiplicative table models with a hydrogen column density of A\va = 5.08 x 10 22 cm -2 
and log ^wa = 0.06. Note that this warm absorber is essentially neutral, implying that it does not 
experience significant heating from the central engine and may exist at some distance from the 
source of the X-ray emission in NGC 4051. Finally, a soft excess is also seen below ~ 1 keV, as 
in so many other Seyfert galaxies. The best fit to this feature is achieved with a thermal blackbody 
component, with kT = 0.27 and a flux ~ 120 times smaller than the hard power-law component. 

A narrow 6.4 keV Gaussian was successfully fit to the core of the Fe-Ka line in this source, 
though by forcing the line to be narrow (i.e., freezing a = 0.0) we note significant residuals mark- 
ing the red and blue wings of the broadened line. The narrow core has an EW = 57.8 eV. The 
remaining broad line is best fit with a kerrconv ( ref 1 ) model having a BH spin of a > 0.67 and 
an r'min = 2.25 — 4.62 r ms . Full parameter values and 90% confidence error bars are listed in Ta- 
ble 10. Fig. 7 shows the iron line residual, best fit to this residual, and relative contributions of the 
individual model components, respectively. 


Aft 7. Ark 120 

Ark 120 is a bright Sy-1 AGN with an estimated BH mass of ~ 2 x 10 8 M© (Wandel et al. 
1999) and a bolometric luminosity of Lb 0 i > 10 45 ergs -1 (Edelson & Malkan 1986). At a red- 
shift of z = 0.0327, this source has a relatively constant 2 — lOkeV X-ray flux of F 2-10 ~ 2.50 x 
10 -1 1 erg cm -2 s -1 (Vaughan et al. 2004). The source is radio-quiet, and due to a lack of observed 
evidence for intrinsic absorption, Ark 120 has been labeled a “bare” Sy-1 nucleus (Ward et al. 
1987; Vaughan et al. 2004). Its host galaxy is an early-type spiral of Hubble type SO/a with an 
inclination angle of i « 26° (Nordgren et al. 1995). Ark 120 has been observed by most of the 
major X-ray observatories. An EXOSAT observation showed the source to have a steep soft X-ray 
spectrum (Turner & Pounds 1989), as did a subsequent ROSAT observation (Brandt et al. 1993). 
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Fig. 7. — Top left: The 2.5 — 10.0 keV spectrum of NGC 4051 fit with a continuum composed 
of two power-law models modified by Galactic photoabsorption, a warm absorber model and a 
blackbody soft excess. Top right: The 2.5 — 10.0 keV best-fit model for NGC 4051, including our 
continuum model and an ionized disk reflection spectrum convolved with our kerrconv relativistic 
smearing kernel. Bottom: A vF v plot depicting the relative flux in each of the model components 
in the spectrum. 
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Model Component 

Parameter 
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A -^-0.64 
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1 
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Table 10: Best-fitting model parameters for the 2.5 — lO.OkeV spectrum of NGC 4051, including 
components and parameter values from 0.6 — 1.5 keV for completeness. The energies from 1.5 — 
2.5 keV were not included due to the presence of absorption edges from the XMM-Newton mirrors. 
Error bars are quoted at 90% confidence. We required the 6.4 keV zgauss line to be intrinsically 
narrow, i.e., o = 0.0. Redshifts were frozen at the cosmological value for the source, in this case 
z = 0.0023. A warm absorber and blackbody soft excess are present in this object, as in many other 
Sy-1 sources. 
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Furthermore, as mentioned above, these X-ray observations showed no indication of any warm 
absorption features. Similar findings were seen in ultraviolet observations (Crenshaw et al. 1999; 
Crenshaw & Kraemer 2001). 

The source was observed by XMM-Newton in August 2003 for an effective duration of 100 ks 
(~ 2.74 x 10 6 photons) by Vaughan et al. , whose results on the EPIC-pn spectrum were reported 
the following year. We follow the data reduction and analysis of these authors, though we restrict 
our attention to the EPIC-pn camera data. Vaughan et al. find that the continuum emission is well 
described by a simple photoabsorbed power-law model with no evidence for complex absorption 
intrinsic to the system. This finding is substantiated by a comparison to the spectrum of 3c273 in its 
radio-loud state, when the source is dominated by jet emission and displays a very flat, featureless 
X-ray spectrum lacking intrinsic emission and absorption features. Vaughan et al. calculated the 
ratio between the spectra of Ark 120 and 3c273, then normalized it by a spectral model for 3c273. 
The resulting spectrum clearly shows that Ark 120, too, possesses a flat, featureless continuum 
above ~ 1 keV and that its only significant emission line is that of Fe-Ka at 6.4 keV. The authors 
use a multiple blackbody component for the soft excess to reduce residuals left over from the 
power-law fit below ~ 1 keV, and a slight curvature of the continuum is noted that is thought to 
indicate the presence of disk reflection in this source (Vaughan et al. 2004). 

Our continuum fit differs from that of Vaughan et al. , most notably over whether or not a 
warm absorber is present at low energies. These authors relied heavily upon the RGS data to 
draw their conclusions about the X-ray spectrum of Ark 120 below ~ 3 keV, whereas we use only 
the EPIC-pn spectrum to produce our global fits to avoid dealing with the calibration problems 
that exist between the two instruments, as discussed in previous Sections. We do find a similar 
photoabsorbed power-law, however: T = 2.25 and N\\ = 1.26 x 10 21 cm -2 (the Galactic value, 
with no evidence for intrinsic cold absorption in the AGN). We detect the presence of a two- 
zone warm absorber in the system with hydrogen column densities of N\ = 1.51 x 10 23 cm~ 2 , 
N2 = 2.90 x 10 22 cm -2 and ionization parameters of log£,i = 2.97 and log £2 = 1-24. Our best 
fit to the soft excess required one blackbody component of kT = 0.13keV and a flux of Fbb = 
1.95 x 10 -6 phcm -2 s -1 , a factor of 60 below the flux of the power-law component. It should be 
noted that adding in the first WA table resulted in a A% 2 /dof = —838 / + 2, and adding in the second 
table further improved the fit by A% 2 /dof = —1163/ + 2. The blackbody component resulted in 
Ax 2 /dof = —4416/ + 2. All of these additions to the global fit are therefore statistically warranted. 

Vaughan et al. detected the presence of broad and narrow components to the 6.4 keV Fe-Ka 
line (EWbroad ~ 100eV, FWnarrow ~ 40 eV). The narrow component could easily be fit with a 
6.4 keV Gaussian profile, as was our Fe-Ka narrow core, but Vaughan et al. found that a diskline 
model worked best for their broad component (r m j„ « 144 r g ). Interestingly, this broad component 
fit to a rest-frame energy of 6.56 keV rather than 6.4 keV, with the latter yielding a substantially 
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worse global fit (Ax 2 = +21.5) (Vaughan et al. 2004). This energy corresponds to the line emission 
being dominated by mildly ionized iron (Fexx-xxn), meaning that resonant trapping and the 
Auger effect should destroy the line (Ross et al. 1996). 

In an effort to solve this puzzle we have fit our iron line profile with two narrow Gaussian 
components, one at 6.4 keV and one at 6.97 keV, representing neutral and highly ionized iron. 
The lines have isWcoid = 54.4 eV and £ Unionized = 19.9 eV. This fit leaves a residual broad feature 
nicely centered at 6.4 keV, which renders this scenario physically consistent with other systems 
and removes the need to explain the presence of the broad line of iron in an intermediate state of 
ionization. We then achieve a best fit to the broad component using a kerrconv(refl) model, 
which also nicely accounts for the additional curvature noted in the spectrum by Vaughan et al. 
From this fit, we have constrained the spin of the BH to be a = 0.63 — 0.68 and rmin = 1.00 — 
1.09r ms . The best-fit parameter values and error bars are presented in Table 11. Fig. 8 shows 
the iron line residual, best fit to this residual, and relative contributions of the individual model 
components for the best fit. 

Though this result is intriguing for its intermediate BH spin value, one must interpret it with 
caution. This source, in particular, seems to have a weak emission feature present in the spectrum 
just below the neutral iron line in energy that sticks out as a “bump” in the spectrum. It is un- 
clear what this feature may represent, and left unmodeled, the spectral fit seems to smooth over it 
rather than fit it properly. This may affect the parameter values and error bars we obtain for the 
kerrconv(refl ) model. More detailed, longer observations with higher spectral resolution will 
be required to solve this puzzle and more reliably determine the BH spin in Ark 120. 


A. 0.8. Fair all 9 

Fairall 9 is a radio-quiet Sy-1 type galaxy with an elliptical companion, both at a moderate 
redshift of z — 0.047. This source has not been observed to undergo very large changes in X- 
ray flux over the time it has been observed: the typical 2 — lOkeV flux of Fairall 9 is F 2-10 ~ 
1.5 — 5.0x 10 -11 erg cm -2 s -1 (Reynolds 1997; Gondoinetal. 2001). Previous X-ray observations 
with ASCA show that the source has a continuum well-described by a photoabsorbed power-law 
with a fluorescent Fe-Ka line and a high-energy tail (a common signature of disk reflection, as 
has been discussed previously with respect to a number of other sources in this Section) (Reynolds 
1997; Nandra et al. 1997). A soft excess has also been detected below ~ 2 keV (Pounds et al. 1994), 
though no strong evidence of warm absorption has been observed in the soft spectrum (Gondoin 
et al. 2001). 

The source was observed with XMM-Newton by Jansen et al. in July 2000 for an effective du- 
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Fig. 8. — Top left: The 2.5 — 10.0 keV spectrum of Ark 120 fit with a continuum composed of a 
power-law model modified by Galactic photoabsorption, a two-zone warm absorber model and a 
blackbody soft excess. Top right: The 2.5 — 10.0 keV best-fit model for Ark 120, including our 
continuum model and an ionized disk reflection spectrum convolved with our kerrconv relativistic 
smearing kernel. Bottom: A vF v plot depicting the relative flux in each of the model components 
in the spectrum. 
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Table 11: Best-fitting model parameters for the 2.5 — 10.0 keV spectrum of Ark 120, including 
components and parameter values from 0.6 — 1.5keV for completeness. The energies from 1.5 — 
2.5 keV were not included due to the presence of absorption edges from the XMM-Newton mirrors. 
Error bars are quoted at 90% confidence. We required the z gauss lines to be intrinsically narrow, 
i.e., c = 0.0. Redshifts were frozen at the cosmological value for the source, in this case z = 0.0327. 
A two- zone warm absorber and blackbody soft excess are present in this object, as in many other 
Sy-1 sources. 
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ration of no 29 ks (~ 3.95 x 10 5 photons) (Jansen et al. 2001), and the EPIC-pn data were presented 
in full by Gondoin et al. the following year (Gondoin et al. 2001). We have followed the general 
procedures for data reduction and analysis discussed by Gondoin et al. , using the most up-to-date 
calibration files and software as we have for all other sources analyzed in this work. 

Gondoin et al. find that continuum of Fairall 9 is best modeled with a photoabsorbed power- 
law (Ah = 3 x 10 20 cm -2 , r = 1.80), in which the absorbing column density of neutral hydrogen 
gas is approximately equal to its Galactic value along the line of sight to this source ( Nn = 3. 19 x 
10 20 cm" 2 ). A soft excess is noted below ~ 2keV, and an RGS analysis yielded a best-fitting 
blackbody model for this component: kT = 0. 17 keV. Incorporating reflection from an ionized disk 
with solar abundances and a reflection fraction of / re fi = 1 .0, the authors find that the high-energy 
tail is quite well modeled by a disk at an inclination angle of i = 26° (Gondoin et al. 2001). We have 
performed a similar fit using only the EPIC-pn data, and find some similar parameters describing 
the power-law, blackbody soft excess and reflection components (the column density of neutral 
hydrogen absorption was frozen at the Galactic value): T = 2.04, kT = 0.20 keV with a flux ~ 0.5% 
that of the power-law component, and the disk inclination angle is i = 57° for solar abundance 
and assumed near-neutrality. Unfortunately, due to the relatively small number of photons, the 
disk ionization and metallicity were not able to be adequately constrained, hence the assumptions 
of solar metallicity and near-neutrality. In contrast to the Gondoin et al. analysis, we do detect 
evidence for a warm absorber within the system with Awa = 1.21 x 10 23 cm -2 and log ^wa = 2.85. 
Inclusion of this component improves the global goodness-of-fit by A% 2 /dof — — 1 3/ + 2, making it 
statistically significant according to the F-test. The presence of the warm absorber also eliminates 
some of the residual continuum curvature in the soft spectrum. 

Gondoin et al.also note the presence of a 6.4 keV Fe-Ka line in the spectrum with EW ~ 
120eV, which they modeled with a Gaussian component. Due to its relatively narrow profile, the 
line is suggested to originate from low-ionization material orbiting relatively far out in the disk. 
Additionally, an Fe-K absorption edge is also noted at 7.64 keV with 1 = 0.18, consistent with 
reflection from cold, optically thick material. We find that the inclusion of such an edge is not 
statistically robust in our model once reflection is included, but we do find the 6.4 keV line in 
emission and model it first with a narrow Gaussian feature representing the core of the Fe-Ka line. 
This feature has an EW = 133 eV, but even after its inclusion residual features remain indicating 
the presence of a broader iron line component. The best fit for this line, which is physically 
consistent with our disk reflection model, is achieved using a kerrconv(refl ) component with 
an inner radius of emission of r m j n < 6.27 r ms , or < 10.75 r g , calculated using the fit value of BH 
spin ( a = 0.972). The BH spin itself, however, was unable to be constrained. Best-fit values and 
error bars for all the parameters in our model are listed in Table 12. Fig. 9 shows the iron line 
residual, best fit to this residual, and relative contributions of the individual model components, 
respectively. 
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Fig. 9. — Top left: The 2.5 — 10.0 keV spectrum of Fairall 9 fit with a continuum composed of a 
power-law model modified by Galactic photoabsorption, a warm absorber model and a blackbody 
soft excess. Top right: The 2.5 — 10.0 keV best-fit model for Fairall 9, including our continuum 
model and an ionized disk reflection spectrum convolved with our kerrconv relativistic smearing 
kernel. Bottom: A vF v plot depicting the relative flux in each of the model components in the 
spectrum. 
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Table 12: Best-fitting model parameters for the 2.5 - 10.0 keV spectrum of Fairall 9, including 
components and parameter values from 0.6 - 1 .5 keV for completeness. The energies from 1 .5 — 
2.5 keV were not included due to the presence of absorption edges from the XMM-Newton mirrors. 
Error bars are quoted at 90% confidence. We required the 6.4 keV zgauss line to be intrinsically 
narrow, i.e., o = 0.0. Redshifts were frozen at the cosmological value for the source, in this case 
z = 0.0470. A warm absorber and blackbody soft excess are present in this object, as in many other 
Sy-1 sources. 
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